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LIBRARY ROUTINE A 1 - 63 



TITLE Floating Decimal Arithmetic Routine 

TYPE Interpretive routine with 18 interpretive orders, entered as 

a closed routine, left by an 8 J interpretive order. 
NUMBER OF WORDS 168 
PURPOSE This routine manipulates numbers in the floating decimal form, 

that is numbers which are represented as A x lO 5 . It is of the interpretive type. 
This means that it selects parameters called interpretive orders which are written 
by the user one at a time and performs a calculation corresponding to each inter- 
pretive order. Interpretive orders carry out normal arithmetic operations such 
as addition and multiplication and some red tape operations such as counting and 
address changing. 

In general, one will use this routine to do computations which 
to not require the full speed of the computer but which are too time consuming to 
be done by hand. It is especially effective for problems with scaling difficulties. 
In a sense one may think of the floating decimal routine as converting the Illiac 
to a medium speed floating decimal computer having a very convenient order code. 

ACCURACY About 9 decimals 

TEMPORARY STORAGE 0, 1, 2 

PRESET PARAMETERS S3 is used to specify two locations of non -temporary storage, 

S3 and 1S3, which are used for the floating decimal accumulator. 

METHOD OF USE The floating decimal routine is entered as a standard subroutine. 

Following the entry, i.e. after the transfer of control to the subroutine, one be- 
gins writing interpretive orders . These orders each occupy one half word and consist 
of a pair of function digits followed by a single address. They therefore have the 
same form as standard machine orders and may be read by the Decimal Order Input 
with full use of the conventional terminating symbols. 

The first of the two function digits of an interpretive order 
describes the group characteristics of the order and may take values 0, 1,..., 8. 
Normal arithmetic interpretive orders have this digit equal to 8. The second of 
the two function digits describes the type of interpretive order. 
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■uramu g fi ra queer list with first fumctioh digit b » 8 

Let F be the floating decimal number la the floating accumulator 
and let F(n) be the floating decinal number In location n. 

80 M Beplace F by F - F(n). 

81 n Replace F by -F(n). 

B2 n Transfer oontrol to the right hand Interpretive order in n if F > 0. 

g«- . Transfer control to the left hand ^jiterpretlve order 5 is n If F ^ 0. 

81* n Beplace F by F + F(n). 

85 n Replace F by F(n). 

86 n Replace F by F/F(n). 

87 n Replace F by F x F(n). 

88 Replace F by one number read from the Input tape punched as sign, 

any number of decimal digits, sign, and two decimal digits to represent 
the exponent* For example, .8971 x 10 would be punched as 

+ 8971 + 10. 

89 n Punch or print F as a sign, n decimal digits, sign, two decimal digits 

to represent the exponent and two spaces . This print out may be re- 
read by this routine* After F has been punched or printed it 
may not remain in the floating accumulator unmodified, n can take 
values 2 to 9* 

8K n Replace F by n if < n < 200 

86 n Replace F(n) by F. 

8H n Replace F by |f| - |F(n)|. 

8j n Transfer control to the ordinary Illlac order on the left hand side of n. 
This used to escape from the floating decimal subroutine. 

6F n Give a carriage return and line feed and start a new block of printing 

having n columns . ' This order is only obeyed once for a particular block 
of printing. At this time a counter is set up which will cause a carriage 
return and line feed to oqcur automatically from then on after every 
set of n numbers that is printed. 
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HBTERPKETIVE ORDERS WITH b ^ 8. If the first function digit of an interpretive 
order is 0, 1^...,7 it will refer to one of a set of control registers or b- 
registers in the floating decimal routine which are similarly numbered. These 
registers are used for counting the number of passages through loops or cycles and 
for advancing addresses on successive passages. For this purpose a particular 
b-registers which may be used in a particular cycle contains two counting indice 
g. and c. . These are both integers in the range to 1023 • The index c is used 
for counting purposes to determine the number of passages through a loop. The 
index g. is used for advancing the addresses of interpretive arithmetic orders. 
Although the interpretive order with first function digit b is not actually altered 
in the memory it is obeyed as if g. were added to its address. The index g, is 
increased by one upon each passage through the cycle. The multiplicity of b-regis- 
ters allows one to program many loops within loops . 

ORDER LIST WITH n ^ 8 

Replace F by F - F(n+g.) 
Replace F by -F(n+g. ) 
Replace g^ c^ by g fe +.1,.:^ + 1. 

Then transfer control to the right hand (if b2 n) or left hand inter- 
pretive (if b3n) order in h if c,+l is negative. This transfer is 
used at the end of a loop. 
Replace F by F + F(n + g,) 
Replace F by F(n + g.) 
Replace F by F/Ffn+g^ 
Replace F by F x F(n+g fe ) 

Replace g., c, by 0, -n. This interpretive order is used for pre- 
paring to cycle around a loop n times. 
Replace F(n + g.) by F 
Replace F by |f| - |F(n+g b )| 

Replace g,, c, by g. + n, c . This interpretive order is used when 
one wishes to step addresses by some increment other than +1 in a 
loop. If one places bL 1022 in a loop the effect will be to decrease 
addresses by one on each passage. bL 1 will increase them by 2 etc. 
8L n Replace g. c, by n, c, , where b is the last b-register referred to by 
some previous interpretive order. 



bO 


n 
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n 
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n 


b5 


n 
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n 


b7 
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bK 


n 
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n 


bH 
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6* 
80 
6k 
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DURATIOH OF IHDIVIDUAL IM ' KKJ tt Ua ' lV E ORDERS 

5 milliseconds + m x (3/2). Where m is the number of shifts re- 
quired to convert A, p back to standard form. 



81 




2. 


— 2 milliseconds 


— 3 milliseconds 


87 


5 milliseconds 


86 


6 millisecons 


8K 


3 milliseconds 


8S 


3 milliseconds 


8F 


3 milliseconds 


8L 


2 milliseconds 


8j 


3 milliseconds 



When an interpretive order is preceeded by b / 8, add one millisecond 
to the above times . 

When one wishes to repeat a cycle of interpretive orders n tines the 
interpretive order bK n may be written before entering the loop to set the counter 
c^ to -i|. The interpretive orders in the loop will be obeyed n times if the loop 
is terminated with a b2 or b3 interpretive order to transfer control to the begin- 
ning of the loop. This transfer of control interpretive order will be obeyed n-1 
times and disobeyed the nth time. 

The following examples illustrate the construction of such loops. 
EXAMPLE 1 Calculate x 10 where x = F(k) 

8K IF Set F » 1 

2K 1QF Set to cycle 10 times 



8? k? c 2 = -10 

23 1L Transfer control to 1 relative 9 times. 



EXAMPLE 2 Replace F(lOO + l) "by F(200 + i) + F(300 + i)^i * 0,1* '",9 
6 



8 



OK 10F 
05 200F 
Ok 500F 
OS 100F 



«o = °' c o - " 10 

F = F(200 + g Q ) 



02 6L 



F = F(200 + g Q ) + F(300 + g Q ) 
F(100 + g )3j F(200 + g Q ) + F(300 + g Q ) 
increase g Q by 1, transfer control 9 times to 
6 relative. 



19 



19-i 



EXAMPLE 3 Evaluate E a ± x y ~ , where a = F(l00 + i), x = F(lO). The operations 



i=0 



form (((ax + a ) x + a ) x + a,) x + -— 



2 


8K F 




OK 20F 


3 


87 10F 




Ok 100F 


k 


03 3L 



Clear F 

Prepare to cycle 20 times 

multiply by x 

add F(100 + g ) 



EXAMPLE k Print 10 numbers each with 5 decimal figures in a block of 3 columns 
from 100 to 109 • 

Start block with 3 columns 
Print 10 numbers 



22 


8F 3F 




OK 10F 


23 


05 100F 




89 5F 


2k 


03 23L 



EXAMPLE 5 



Place the number F(200 + i) in 300 + 2i for i * 0,1,... ,2k. 



5 


OK 25F 




IK F 





05 200F 




IS 300 F 


7 


1L 2F 




03 6l 



g - g-L = 
c Q = -25 



advance g by 2 
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EXAMPLE 6 H*ce the number F(200 + l) in 300 - i for i « 0,1,... ,24 



5 


0K25P 




1KF 


6 


05 20QF 




IS 50QF 


7 


1L 1023F 




03 6L 



reduce g, by 1 
increase ^ and c by 1 

Use of Auxi liary Routines . It is often convenient to be able to 
leave the floating decimal routine so as to modify interpretive orders or to 
perform calculations which may be done more effectively outside of floating point. 
aw jiVsYc »n™ i xoa wing decimal routine one uses an 8J n oT^dsi* fAii n4 , .AT%f4nw4 
floating decimal auxiliaries are entered in this way.) !Po return to floating decimal 
one should transfer control to the left hand side of word 29 of the floating decimal 
routine. The interpretive order following the Bj n order which was last obeved will 
then be obeyed and so on. In this way it is not necessary to plant a link in auxi- 
liary subroutines. One may, in fact, think of the 8j n order as a subroutine order. 
In case any changes are made in the floating decimal accumulator while outside the 
floating decimal routine, control should be returned to the left hand side of word 
19 rather than 29 so that this number may be standardized before reentry. 

Handling o f lumbers Each number is represented in the foimA x 1Q 9 
where 1 > |a| > l/lO, and 64 > p > -64. in a single register of the memory the 
number A is placed in the 33 most significant binary digits (a , a ,,.., a ) 
in the same way as an ordinary fraction is placed in the entire register. An 
accuracy of between 8 and 9 decimal digits is therefore achieved. The exponent 
p is stored as the integer p + 64 in the 7 least significant digits of the same 
register. For convenience the floating decimal accumulator uses two registers 
S3 and 1S3 for holding the number A x 10*. The fraction A/2 is in S3 and the 
integer p + 64 is in 1S3. 

The only exeption to the above rules is the number zero which can- 
not, of course, be represented as A x 10* with |a| > 10. For this reason zero is 
handled in a special way. It is represented as a number with A - and p * -64. 
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This representation happens to correspond exactly with the ordinary machine 
representation of zero. 

After each arithmetic interpretive order is obeyed the number in the 
floating decimal accumulator is standardized, i.e# the number of S3 representing 
A/2 is adjusted so 1 > |A| > l/lO and p is changed accordingly. To accomplish 
this control is transferred to word 19 in the floating decimal routine after each 
arithmetic order. 

If an interpretive store order is attempted when F has an exponent 
greater than 63 the machine will stop on the order 5^ P at location p, where p is 
word 72 of the routine. 

Important Words in the Routine . Word 2 in the floating decimal rou- 
tine determines the location of the current interpretive order. When obeying the 
left hand interpretive order in location n this word is 50 nF S5 20F and when 
obeying the right hand interpretive order in location n it is L5 nF 00 20F. Other 
words of interest are the b-registers which start at word 158 (for g-> and c ) and 
go to 165 (g 7 and c 7 ). These register hold g, and c, in the form 

80 g b F 00 (2048 + c^) F. 

Warning When the same number is continually added to a sum, such 
as when an argument is being increased, the error can be quite large, because it 

is additive over a decade. For example, if we increase 10 to 100 by units we can 

-33 
get a maximum error of 90 x 2 because the errors all have the same sign. If we 

3 k -33 

increase 10 to 10 we can have a maxlTmrm error of 9^000 x 2 ^ . This can easily 

be prevented by writing an auxiliary subroutine to stabilize the fractional part of 

F, i.e. to replace it by the nearest multiple of say 10 . Such a subroutine could 

be as follows 

m 50 S3 

A contains nearest multiple of 10" ' 

Location 1 contains zero 

converts to fraction again 

return to 18th order of Routine A-l 

m + 3 I 00 F 

I 00 20,000,000F 2 x 10 7 x 2" 39 ^ q/^ 



m 


50 S3 




7J m+3F 


m + 1 


50 IF 




66 m+3F 


m + 2 


S5 F 




22 S+18F 



DATE November IB, 1^5 



BY David J Wheeler approved BY J. P. Hash 



LOCATION 
NOTE: 



ORDER 



3 

4 

5 

6 

7 
8 

9 

10 
11 
12 

13 
14 



00 K(A1) 

00 59F 
LO 166L 
LO 124L 
40 2L 
50 (n)F 
S5 20F 
40 2F 
32 150L 
46 6L 
10 12F 
L4 125L 
46 9L 
L5 (m)F 
10 IF 
40 F 

10 6F 

01 7F 
40 IF 
26 F 
40 F 

L5 1S3 
LO IF 
50 F 
32 15L 
40 2F 
S5 F 
50 S3 
40 S3 
L5 IP 
40 1S3 



NOTES 
Assume N(S3) * k/2 where F « A x 10* 

N(1S3) = p + 64 

F(m) * B x 10 q 
new value of quantity is indicated "by a 
dash. 



} 



Special orders for entry 

Plant either 50 nF S5 20F or L5 nF 0020F 

in 2L 



„ Select order to be obeyed 

Adjust address if necessary 
Insert address in 6L 

- Form switch order 



- Select argument 
Place fraction B/2 in F 



and exponent q in IF 



\ 



Test exponent difference p - q 



- Interchange A and B 



and place q in 1S3 



LOCATION 


ORDER 


NOTES page 2 


15 


LI 2F 








L4 125L 






16 


k2 17L 

LO 126L 




Insert exponent difference in 17L 


17 


36 29L 




Skip addition if exponent difference 




7J (31 + |p-q 


|)F 


|p - q| > 10 


18 


Lk S3 




-A + B 10" P+<1 or B + A uf*™ 




40 S3 








19 


LL S3 
36 27L 




is A 1 > 1/2 




20 


50 S3 
7J 32L 




Then multiply A f by l/lO 




21 


^0 S3 
L5 ^3L 




- 


Standardize A 
and P so that 


22 


Lk 1S3 
kO 1S3 




and increase p' "by 1 


1 > |A| > 1/10 


23 


26 27L 
50 S3 








2k 


75 127L 
00 kF 




A' x 10 




25 


kO S3 
LI 43L 








26 

1 

i 


Lk 1S3 ' 
kO IS 3 




~p - 1 




27 


L3 S3 

36 68L 




Test if |A ! | < l/lO 
or equal to zero 




28 


Lk 42L 
33 23L 








29 


L5 2L 




if N(2L) is positive (50 nF S5 20F) 




36 1L 




select right hand order in n. 


30 


Lk 128L 
22 1L 




otherwise select left hand order in n+1 


31 


7L 4095F 








LL 4095F 




10 - almost 



LOCATION 
32 

33 

3^ 

35 

36 

37 

33 

39 

to 

1*2 
1*1* 

hi 
1*8 



ORDER 

OOF 00 1000 
0000 0000 J 
OOF 00 100 
0000 0003 J 
OOF 00 10 
0000 0000 J 
OOF 00 1 

0000 0000 J 
OOF 00 
1000 0000 J 
OOF 00 

100 0000 J 
OOF 00 
10 0000 J 
OOF 00 

1 0000 J 
OOF 00 
1000 J 
OOF 

00 61*F 
OOF 00 1*99 
9999 9998 J 
00 F 
00 IF 
LI F 
22 9L 
LI F 
26 67L 
L5 131L 
26 133L 
L5 131L 
22 133L 
JO 1021F 
26 10L 



BOTES 
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10 



-1 



10 



-2 



10 



-3 



10 



-1* 



10 



-5 



10 



-6 



10 



-7 



10 



-8 



10 



-9 



} 



bO 

bl 
b2 



^3 

4)1* also used as a 
binary switch' 



Table of switch 
orders to transfer 
control to the 
appropriate sets of 
orders • 



LOCATION 


ORDER 


IOTES 


mGE k 


*9 


L5 F 
26 67L 




b5 




50 


LI IF 
26 63L 




b6 




51 


LO klL 
26 60L 




b7 




52 


k9 IF 
22 77L 




b8 




53 


L5 6l 
22 97L 




b9 


> 


5^ 


kl IF 
26 138L 




bK 




55 


».L5 6L 
26 69L 




bS 




56 


LT S3 
26 76L 




bN 




57 


L5 2F 
22 5L 




bj 




58 


92 131F 
26 95L 




bF 




■i 59 


L5 6L 
22 lk9L 




bL 




60 


1^ 1S3 
kO 1S3. 




p' »p + q 




6l 


50 F 
7J S3 




A'»AB 


m Multiplication 


62 


00 IF 








63 


22 18L 
IX 130L 
IM 1S3 




p» » p . q + 1 




6k 


40 1S3 
50 42L 






^ Division 


65 


75 S3 




^ 






66 F 




A '' Ax (VlO) + B 





LOCATION 


OHDER 
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66 


35 31L 






61 


22 18L 
kO S3 




p 1 = q 






L5 IF 




_ Clear add* 

_ 


68 


kO 1S3 
26 2% 




A' - -B 


Subtract 


69 


k6 75L 






70 


32 71L 








_ 1 
Test if p + 6k is negative 






kl S3 








71 


kl 1S3 
LO 167L 




„ „ 




72 


3k 72L 
L5 S3 




Stop machine if p > 6k. 




73 


10 6f 
50 IF 




Form packed A, p 


■Store order 


7* 


00 7F 








75 


IX 1S3 
kO (n)F 




and place in m. 










26 29L 






76 


40 S3 




A' = |A| 




L3 F 




B* * -|B| Modulus order 


77 


22 % 
81 kF 




t 


78 


00 39F 








36 80L 




Input Sign place + l/2 in IF 


79 


LI IF 








kO IF 






80 


27 8^L 
50 IF 






81 


7J 32L 
kO IF 






82 


50 IF 
75 2F 







LOCATION 
83 
84 

85 
86 

87 
88 

89 
90 
91 
92 
93 
9^ 
95 
96 

97 
98 

99 



ORDER 

00 39F 
L4 S3 
40 S3 
81 4f 
40 2F 
LO 129L 
36 87L 
22 80L 
81 4F 
40 F 
8l 4F 

50 F 

74 129L 
LI 2F 
L4 129L 
36 92L 

51 31L 
22 92L 
S5 F 
L4 4lL 
40 1S3 
26 19L 
LI ( )F 
40 4f 
L5 2F 
10 20F 
42 94L 
47 94L 
26 29L 
L4 91L 
46 105L 
L5 94L 
LO 79L 
36 101L 



NOTES 



Input conversion cycle 
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Read 



Exponent 



T- Control to 19L to standardize A and p 



Column Count 

Arrange new column count and 
give a carriage return and 
line feed. 

Set address of round-off order 



Give a carriage return 



LOCATION 


ORDER 
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100 


92 131F 
22 122L 




When necessary 




101 


k6 9kL 
L5 S3 








102 


Ilq p 
36 104L 






Print order 


103 


92 708F 
22 lO^L 








10^ 


92 6kkF 
50 101L 




-Print sign and round -off 


■ 


105 


7J ( )F 
L6 F 








106 ■ 


l6 f 








32 110L 




Test if round-off + |a| > 1 


107 


L5 ^3L 
IA 1S3 






108 


ifO 1S3 




If so replace A by l/20 




L5 ^2L 




— 


109 


kO S3 




and p by p + 1 




kO F 






110 


22 lO^L 








10 39F 






111 


75 129L 








00 36F 




Print cycle 


112 


82 kF 








10 1+OF 






113 


L5 6l 
LO k9L 






114 


k6 6l 
LO 62L 






115 


36 111L 








92 96lF 




Space 


116 


LI kQL 




Binary switch 


/ 


kO k&L 




- ■— ■ - 



LOCATION 


ORDER 
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117 


32 118L 




Did we print fraction or exponent? 




92 96IF 




Space 


118 


26 29L 








k6 6l 






119 


L5 IS 3 
LO to 




Arrange 00 prxn^ exponenw 


120 


kO F 








50 F 




- 


121 


75 33L 
00 38F 






122 


26 102L 








92 513F 






123 


00 2^F 
26 101L 






12k 


5S F 




L5 (n)F 00 20F - 50 (n)F S5 20F 




S5 F 






125 


00 172L 




" Base of switch order = kk + 128 




00 31L 




Base of lO* 5 table order 


126 


00 172L 




""[End test for 10" p table order 




00 42L 






127 


50 F 




_10/l6 




00 F 






128 


5S IF 








S5 F 




_ 50 (n+l)F S5 20F = L5 (n)F 00 20F 


129 


00 F 
00 10F 




-39 
10 x 2 Jy 


130 


00 F 

00 65F 




65 x 2" 59 


131 


50 F 




^Base for 2,3 orders 




S5 20F 






132 


00 IF 




_Double unit for counting 




00 IF 






133 


LO 12^L 










kO IF 




-Prepare for transfer 


. 



LOCATION 


ORDER 
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134 


L5 2F 
46 IF 








135 


32 l43L 






_ orders 




L5 S3 




Transfer control if positive 


b3 


136 


36 137L 
26 2% 




Select next order if negative 




137 


id IF 
22 IL 






138 


L5 2F 




Is b * 8 or £ 8? 




36 146L 








139 


46 IF 
50 IF 




Place m/200 in S3 




l4o 


75 33L 
00 loF 






8K 


l4l 


40 S3 
L5 130L 




and 66 x 2* 59 in 1S3 




142 


L4 43L 




■ ' 




40 1S3 






*43 


26 19L 




■ 






L5 (157 + b)L by 152 


b J 8 


144 


L4 132L I 


Increase c. and g. by 1 


r b2. . 

orders 




40 (I57+b)L 




b3 


145 


00 28F 
26 I36L 




Then test sign 


146 


10 2QF 
42 IF 






bKb + 8 


Ikf 


U IF 
50 156X 






_make g. = 
c, = 2048-n 


148- 


L4 157L 










40 .,{157+1 


»)L 






149 


26 2% 








46 (l57+b)L 


From OL switch 


150 


26 29L 








50 2F 







LOCATION 


ORDER 
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151 


01 3F 
L4 ihjL 






152 


42 l43L 
42 155L 




Adjust all addresses 


153 


42 148L 
42 149L 




dependent on b 


15^ 


42 144L 
JO 166L 






155 


11 3F 
L5 157L 




Modify order add g, + 1 


156 


S4 F 
26 4L 






157 


80 F 

00 2048F 




Constants 


158 


80 F 
00 F 






159 


80 F 
OOF 






160 


80 F 
00 F 






161 


80 F 
00 F 






162 


80 F 
00 F 






163 


80 F 
00 F 






164 


80 F 

00 F 






. 165 


80 F 
00 F 






166 


LL 4094F 
4K 4076F 






167 


00 F 
00 128F 







M 



